Experimental studies of equilibrium vortex properties in a Bose-condensed gas 
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We characterize several equilibrium vortex effects in a rotating Bose-Einstein condensate. Specif- 
ically we attempt precision measurements of vortex lattice spacing and the vortex core size over 
a range of condensate densities and rotation rates. These measurements are supplemented by nu- 
merical simulations, and both experimental and numerical data are compared to theory. Finally, 
we study the effect of the centrifugal weakening of the trapping spring constants on the critical 
temperature for quantum degeneracy and the effects of finite temperature on vortex contrast. 
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I. INTRODUCTION 

After the initial observations of vortex lattices in Bose- 
Einstein condensed gases Q, 0, 0, 0] , most of the exper- 
imental work has focused on dynamical behavior of vor- 
tices and lattices, including Kelvons HUSH; Tkachenko 
waves d, H EH jljl [3 HI, and various nonequilib- 
rium effects |l4l llfjj . Equilibrium properties, in con- 
trast, have been relatively neglected by experimenters. 
This imbalance is not indicative of a lack of interesting 
physics in equilibrium behavior, but simply reflects the 
usual experimentalist's preference for measuring spectra 
rather than static structure. Theorists, on the other 
hand, have investigated equilibrium prop erties exten- 

siveiy [IlIIlIillilliaBlillilliilllllllli, 

and our purpose in this paper is to partially redress this 
imbalance with a series of experimental studies focusing 
on equilibrium properties of rotating condensates. 

The vortex lattice in a rotating Bose-condensed gas 
naturally organizes into a regular triangular lattice, or 
Abrikosov lattice, originally observed in superconduc- 
tors. The lattice can be well characterized by the nearest- 
neighbor lattice spacing and by the radius of each vortex 
core (b and r v respectively for the purpose of this pa- 
per). The nearest-neighbor lattice spacing, 6, is generally 
thought to be determined only by the rotation rate when 
in the high-rotation regime where the rotating BEC ex- 
hibits nearly rigid-body behavior. Numerical work 
and early analytical work however, suggests that this 
rigid-body assumption yields lattice constants that are 
smaller than would be seen in the case of a finite-size 
trap ped BEC. Recent work by Sheehy and Radzihovsky 
|l7| has tackled analytically this discrepancy and found 
it to be a necessary consequence of the inhomogeneous 
density profile of the condensate. With this theory they 
address the question of why the lattice is so remark- 
ably regular given the condensate density profile. They 
also derive a small, position-dependent, inhomogeneity- 
induced correction term to the lattice spacing. An in- 
teresting implication of this theory is that the vortices 
must move slightly faster than the surrounding superfluid 



even near the rigid-body limit. More striking still is the 
prediction that the superfluid should exhibit a radially- 
dependant angular velocity (or radial shear flow), which 
directly follows from their calculation of inhomogeneous 
vortex density. While a differential rotation rate is not 
directly observable in our system, the position-dependent 
variation of the nearest-neighbor lattice spacing is stud- 
ied in Sec. II below. It should also be noted that the 
inhomogeneity in the areal density of vortices, predicted 
in Ref . [l2| , can also be derived in the limit of the lowest 
Landau level (LLL). This property of the LLL was first 
brought to our attention by A.H. MacDonald and has 
been the subject of two recent publications by Watanabe 
et al. and Cooper et al. 

The second effect we study in this paper concerns the 
core size of the vortices. Once rotation rate and den- 
sity are fixed, the vortex core size is a length scale that 
the condensate chooses on its own. In this sense vor- 
tex core size constitutes a fundamental property of the 
system and has therefore been the subject of much the- 
oretical work By analogy to superfluid 4 He, 
the core size is dictated by the atomic interactions and 
is of order of the healing length. For our system the 
healing length is only one and a half times the average 
interatomic spacing. Because of this diluteness one might 
wonder if there are certain regimes of sufficiently low or 
high density where one would see a deviation from mean- 
field theory. Investigation of core size makes up Sec. Ill 
of this paper. 

It has been predicted 23] that at higher rotation the 
condensate begins to enter the LLL regime where the 
condensate wave function is constrained to occupy only 
the LLL harmonic oscillator states. Our condensate is 
strikingly different from quantum Hall systems also asso- 
ciated with the LLL in that the mean-field approximation 
is still a valid way to account for interactions. In this case 
the vortex core-size to lattice-spacing ratio saturates at 
the LLL limit IH El As Baym and Pethick ^ note 
this requires the cores to deviate from the Thomas- Fermi 
prediction in shape and size. 

Employing a numerically generated Gross-Pitaevskii 
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(GP) wave function, we study this transition and monitor 
the core size and shape over a range of conditions. These 
simulations show a smooth transition from the Thomas- 
Fermi regime to the LLL regime. Additionally we use 
numerical analysis to study the shape of the vortex cores 
and verify that when the core-area saturation occurs the 
vortex-core wave function is the one predicted for the 
LLL regime. This result suggests that the fractional core 
area is a possible means to probe the transition to the 
LLL regime. To test this idea, we compare simulations 
with experimental data. This comparison and possible 
systematic errors in these experimental measurements 
are discussed in Sec. Ill B. 

In Sec. IV we examine the rotational suppression of 
quantum degeneracy. This effect is due to centrifu- 
gal forces weakening the radial trap-spring constants. 
Weaker spring constants then lead, in a straight forward 
way, to a lower critical temperature for fixed numbers of 
atoms. 

Finally in Sec. V we examine a proposal that a mea- 
surement of the contrast of vortex cores could serve as a 
sensitive thermometer for a condensate in the regime for 
which the temperature is less than the chemical potential 
and other methods of thermometry become unreliable. In 
Sec. V we discuss our preliminary efforts to realize this 
vision. We are able to see an effect, but we have not yet 
been able to extend this measurement technique below 
the usual limits. 

The rest of Sec. I discusses experimental issues com- 
mon to all the results of this paper 

A. Experiment 

Our experiment begins with a magnetically trapped 
cloud containing greater than 10 7 87 Rb atoms in the 
\F = l,mp = —1) hyperfine ground state, cooled to 
a temperature roughly three times the critical tempera- 
ture (T c ) for Bose condensation. Using a TOP trap we 
confine these atoms in an axially symmetric, oblate and 
harmonic potential with trapping frequencies ({u) p , lu z } — 
27r{7, 13}Hz) with axis of symmetry along the vertical (z) 
axis. Rotation is generated in the thermal cloud by res- 
onantly coupling to a scissors mode of the cloud [Hi ED • 
To do this we gradually apply an elliptical deformation 
to the magnetic trapping potential by distorting the am- 
plitude of the rotating TOP field in time. The resulting 
distorted trap has roughly similar average radial-trapping 
frequency but an ellipticity 32] in the horizontal plane 
of 25%. The uncondensed cloud is held in this trap for 
5 seconds while any excitations die out. At this point 
the angle of the major axis of the elliptically distorted 
trapping field is jumped quickly by 45 degrees in the 
horizontal plane to generate the initial conditions of the 
scissors mode. From here the cloud is allowed to evolve 
for 155 ms, or roughly one quarter period of the result- 



ing scissors mode, at which point we transfer the cloud 
to a radially symmetric trap. Essentially we have caught 
the scissors oscillation between turning points where all 
the initial linear velocity has turned into rotational veloc- 
ity. Using this method we can generate a cloud rotating 
at roughly half the radial trap frequency, with minimal 
heating. By lowering the amplitude of the trap distortion 
or the angle by which the trap is jumped, we can easily 
generate more slowly rotating clouds as well. 

At this point we begin a second phase of rf evapora- 
tion, but this time we evaporate in one dimension along 
the axis of rotation Q. The motivation for this seem- 
ingly inefficient technique is that the ID evaporation al- 
lows us to remove energy from the z axis of the conden- 
sate without preferentially removing high angular mo- 
mentum atoms. Lowering the energy per particle with- 
out lowering the angular momentum per particle accel- 
erates the cloud rotation rate Q. To perform this ID 
evaporation, we adiabatically ramp to a prolate geom- 
etry ({u! p ,u! z } = 27r{8.3, 5.3}Hz) where the rest of the 
experiment is carried out. 

Reaching significant rotation rates by the end of evap- 
oration requires that the lifetime of the thermal cloud's 
angular momentum be comparable to the evaporation 
time. The nearly one-dimensional nature of the evap- 
oration together with the low average trap frequencies 
makes cooling to BEC in the prolate trap very slow (2 
minutes). We obtain angular momentum lifetimes this 
long by shimming the TOP trap's rotating bias field to 
cancel the few percent azimuthal trap asymmetry that 
exists despite careful construction. With this technique 
we suppress the azimuthal trap ellipticity to less than one 
part in a thousand. 

After the evaporation we have a condensate with 
as many as 4.5 million atoms and rotation rates from 
fl = (0 — 0.975)w p , with no observable thermal cloud. 
Rotation can be accurately determined by comparing 
the condensate aspect ratio A (defined as the axial 
Thomas-Fermi radius over the radial Thomas-Fermi ra- 
dius R z /Rp) to the trap aspect ratio Ao = (lu p /lu z ), and 
using the now standard relation 

n/w p = V(i-(a/a ) 2 ) . (1) 

B. Expansion 

The clouds described here typically contain between 
1-200 vortices, each of which is too small to be observed 
in-trap but can be seen after expansion of the cloud. Our 
expansion technique, while not completely unique, is un- 
usual enough to warrant a description. For the exper- 
iments presented here and elsewhere, we need a large 
radial expansion to make sure that the vortex cores are 
large compared to our imaging resolution. Additionally 
we need to suppress the axial expansion in order to pre- 
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serve certain length scales in the condensate as discussed 
below. Clearly, given our low and nearly isotropic trap 
frequencies, the usual expansion technique of shutting 
off the magnetic field and dropping the cloud would not 
meet such requirements. The solution, which has been 
demonstrated by other groups |3^| , is to perform an anti- 
trapped expansion. Rather than simply shutting off the 
trapping potential, we invert the trap in the radial di- 
rection so that the cloud is actively pulled apart. Si- 
multaneously, the magnetic field gradient in the vertical 
direction is used to support against gravity. Because it 
has generated some interest we describe this technique 
in the following paragraph in what otherwise might be 
excessive detail. 

The expansion is achieved in several steps that take 
place in rapid succession. First we employ a microwave 
adiabatic rapid passage technique (ARP) to transfer the 
atoms from the weak-field seeking \F = l,mp = — 1) 
state to the strong-field seeking \F = 2,mp — —1) state. 
The microwave field employed is powerful enough to per- 
form the transfer in lOfis but, as will be discussed shortly, 
we often take as long as 300/iS for this transfer. After 
transfer to the anti-trapped state, the cloud still sits in its 
original position below the quadrupole zero, which means 
that both gravity and the magnetic field are acting to pull 
it downward. To counter this force a downward uniform 
vertical magnetic field is added to pull the quadrupole 
zero below the condensate so that the magnetic field gra- 
dient again cancels gravity. The field is applied within 
10/iS, fast compared to relevant time scales. In this man- 
ner the cloud is again supported against gravity. To re- 
duce curvature in the z direction, the TOP trap's rotating 
bias field is turned off leaving only the linear magnetic 
gradient of the quadrupole field. This gradient is tuned 
slightly to cancel gravity. 

Using this technique, we are able to radially expand the 
cloud by more than a factor of 10 while, at the same time, 
seeing less than a factor of two axial expansion. Unfortu- 
nately even this much axial expansion is unacceptable in 
some cases. In the limit of adiabatic expansion, this fac- 
tor of two decrease in condensate density would lead to 
an additional \[2 increase in healing length during expan- 
sion. Thus, features that scale with healing length, such 
as vortex core radius in the slow rotation limit, would be- 
come distorted. The effect of axial expansion on vortex 
size was first noted by Dalfovo and Modugno [34| . 

To suppress the axial expansion, we give the conden- 
sate an initial inward or compressional impulse along the 
axial direction. This is done by slowing down the rate 
at which we transfer the atoms into the anti-trapped 
state. The configuration of the ARP is such that it trans- 
fers atoms at the top of the cloud first and moves down 
through the cloud at a linear rate. These upper atoms 
are then pulled downward with a force of 2g (gravity 
plus magnetic potential), thus giving them an initial in- 
ward impulse. Finally, the ARP sweep passes resonantly 



through the lowest atoms in the cloud: they, too, feel a 
downward acceleration but the axial magnetic field gra- 
dient is reversed before they can accumulate much down- 
ward velocity. On average the cloud experiences a down- 
ward impulse, but also an axial inward impulse. Nor- 
mally the ARP happens much too fast for the effect to 
be observable but when the transfer time is slowed to 
200 — 300^ts the effect is enough to cause the cloud to 
compress axially by 10-40% for the first quarter of the 
radial expansion duration. The cloud then expands back 
to its original axial size by the end of the radial expan- 
sion. 

Despite our best efforts to null out axial expansion, we 
observe that the cloud experiences somewhere between 
20% axial compression to 20% axial expansion at the time 
of the image, which should be, at most, a 10% systematic 
error on measured vortex core radius. The overall effect 
of axial expansion can be seen in Fig.^ where image (b) 
and (c) are similar condensates and differ primarily in 
that (c) has undergone a factor of 3 in axial expansion 
while in (b) axial expansion has been suppressed. The 
effect on the vortex core size is clearly visible. 

Because almost all the data presented in this paper 
is extracted from images acquired after the condensate 
expands, it is worth discussing the effect of radial ex- 
pansion on the density structure in the cloud. In the 
Thomas-Fermi limit, it is easy to show that the anti- 
trapped expansion in a parabolic trap, combined with 
the mean-field and centrifugally driven expansion of the 
rotating cloud, leads to a simple scaling of the linear size 
|3fil | of the smoothed, inverted-parabolic density envelope. 
As R p increases, what happens to the vortex-core size? 
There are two limits that are easy to understand. In 
a purely 2D expansion (in which the axial size remains 
constant), the density at any spot in the condensate co- 
moving with the expansion goes as 1/-Rp, and the local 
healing length £ then increases over time linearly with 
the increase in R p . In equilibrium, the vortex core size 
scales linearly with £. The time scale for the vortex core 
size to adjust is given by hj /i where fj, is the chemical po- 
tential. In the limit (which holds early in the expansion 
process) where the fractional change in R p is small in a 
time h/ /j,, the vortex core can adiabatically adjust to the 
increase in £, and the ratio of the core-size to R p should 
remain fixed as the cloud expands. 

In the opposite limit, which applies when R p expands 
very rapidly compared to H//J,, the inverted parabolic po- 
tential dominates the dynamics, and every point in the 
cloud expands radially outward at a rate proportional to 
its distance from the cloud center. In this limit, the "fab- 
ric of the universe" is simply stretched outward, and all 
density features, including vortex core size expand at the 
same fractional rate. Again, the ratio of core-size to R p 
should remain fixed. 

So in the two extreme limits, the ration of vortex core 
size to R p (and other density features, such as nearest- 
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FIG. 1: Examples of the condensates used in the experiment 
viewed after expansion. Image (a) is a slowly rotating con- 
densate. Images (b) and (c) are of rapidly-rotating conden- 
sates with similar in-trap conditions. They differ only in that 
(c) was allowed to expand axially during the anti-trapped 
expansion. The effect on the vortex core size is visible by 
eye. Images (a) and (b) are of the regularity required for the 
nearest-neighbor lattice spacing measurements. 



neighbor vortex separation) remains fixed. It is reason- 
able then to assume that this behavior will be true in 
general in the intermediate regime between the two ex- 
pansion rates. Extensive numerical simulations were per- 
formed to validate this assumption. 

On a separate note, the vortex lattice spacing, unlike 
vortex core size, is largely unaffected by the presence of 
axial expansion. However, we see some indication that 
deep in the lowest Landau level large axial expansions 
can affect vortex spacing as well as size. 

Once expanded, the cloud is imaged along the vertical 
direction, and data is extracted by fitting the integrated 
(along the line of sight) condensate density profile with 
a Thomas-Fermi distribution. We then subtract this fit 
from the image and easily fit the remaining vortex-core 
residuals with individual 2D Gaussian profiles to deter- 
mine the core centers and radii. For the purpose of this 
paper the vortex radius r v is defined to be the RMS ra- 
dius of the 2D Gaussian that we fit to the core. For 
clarification, the RMS radius of a 2D Gaussian would be 
0.60 times its FWHM. Condensate and vortex fits can be 
performed iteratively to reduce error. 

Before each expanded image we also take a horizon- 
tal, nondestructive, in-trap image of the cloud immedi- 
ately before expansion. From this image, rotation rate 
and atom number are determined. Length scales in the 
expanded cloud can be scaled back to in-trap values 
by dividing by the radial expansion factor, defined as 
R p (expanded) / R p (in — trap). 



C. Numerical studies 

The numerical studies discussed in this paper are done 
by setting up an initial in-trap condensate wave function 
with a given N and f2 on a 2048x2048 lattice and then 
relaxing this wave function by propagating the Gross- 
Pitaevskii equation in imaginary time. All work shown 
in Sec. Ill B and III C is done in 2D. Additionally a 
radially symmetric 3D simulation can be performed for 
a single vortex, as is done in Sec. Ill A. Once the final 
wave function is found, we convert to an atom density 
profile which can be fit and analyzed in the same manner 
as the experimental data. 



II. THE LATTICE CONSTANT 

At first sight, vortex lattices, such as the one seen in 
Fig. ^b), appear perfectly regular. However as noted in 
the introduction, Sheehy and Radzihovsky ^3] predict 
that there should exist a small correction to the vortex 
density in the condensate due to the condensate, density 
inhomogeneity. One result from Ref. 01 is thai the 
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areal density of vortices is 



n v (p) 



1 



2kri{\ - ( P /R p y 



T ln[h/ '(2.718 mnS, 2 )} 



(2) 

where m is the mass of rubidium and £ is the healing 
length (calculated from the measured density). This 
equation can conveniently be thought of as the rigid body 
rotation (first term) plus the density inhomogeneity cor- 
rection that reduces vortex density. We compare to ex- 
perimental measurements by converting vortex density 
to a nearest-neighbor lattice spacing, conveniently ex- 
pressed in units of condensate radius 



b(p) = ^2/(31/2^))-! 



(3) 



To study this lattice inhomogeneity effect experimen- 
tally, we generate condensates with rotation rates of 
Q/ujp = 0.5 — 0.9. To extract the vortex separation, we 
expand the cloud by a factor of 10 in the radial direction 
using the anti-trapped expansion technique. The con- 
densate and vortices are fit as described in Sec. I B. The 
nearest-neighbor separation for a given vortex is mea- 
sured by averaging the distance from the vortex center 
to the centers of the six nearest vortices. Because of low 
signal, vortices further than 0.9 R p from the condensate 
center are disregarded. Any remaining vortex with fewer 
than six nearest neighbors (i.e., a vortex in the outer 
ring) is used as a neighbor to other vortices but is not 
itself included in the final data. Obviously using the six 
nearest neighbors assumes a triangular lattice structure, 
so before fitting, each image is checked for defects in the 
lattice. Any image exhibiting broken lattice planes is 
not considered. Once the nearest-neighbor separation is 
measured, it is normalized by the expanded condensate 
radius to compare to Eq. |3| For this comparison, R p , f2, 
and £ are measured or calculated from an in-trap image. 
To improve the theory fit, we allowed R p to be an ad- 
justable theory parameter, but, in each case, the theory 
fit value for R p was within 5% of the experimentally de- 
termined value of R p . This < 5% difference is within the 
calibration uncertainty for such a measurement. Noise is 
suppressed by taking an azimuthal average of the lattice 
spacing data. Due to the discrete nature of the vortices 
this is equivalent to binning the lattice-spacing data by 
radial displacement of the vortex from the condensate 
center. 

Figure [3 shows a comparison to theory for three phys- 
ical condensates and one numerically generated conden- 
sate density profile. Figure Ufa) is data taken from the 
condensate in Fig. Ofa). The two points shown corre- 
spond to the measured vortex density for the center vor- 
tex (first point) and the average vortex density for the 
first ring of vortices. Also plotted is Eq.|3 (solid line) and 
the expected nearest-neighbor lattice spacing for rigid 
body rotation (dashed line). The imperfect fit may be 
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FIG. 2: Measured and binned lattice spacings as a function of 
radial position p. The solid curve is the theory result (Eq. ©) 
of Sheehy and Radzihovsky [l7|. The rigid-body-rotation 
rate lattice spacing is also plotted for comparison (dashed 
line) . Plots (a-c) are experimental data with increasing rota- 
tion. Plot (a) and (c) are data taken from the condensate in 
Figs. 0^a) and (b), respectively. Plot (d) is the same effect 
observed in the numerical data. One can see that theory and 
experiment show a similar dependence on radial position and 
that the fractional amplitude of the density inhomogeneity 
effect is suppressed at higher rotation. Plot (e) is the data in 
(c) plotted without suppressing the zero. The vortex lattice 
spacing changes less than 2% over a region in which the atom 
density varies by 35%. 
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partly due to the discrete nature of the data, vis-a-vis a 
continuum theory . Plots (b) and (c) are condensates 
with increasing rotation rates where (c) is taken from im- 
age mb). Plot (d) is a comparison with numerical data 
prepared with parameters similar to the experimental sit- 
uation in (c). Figure [2fe) is the same data in Fig. ^c) 
but plotted without the suppressed zero to emphasize 
the smallness of the position dependant effect. The areal 
density of vortices is constant to 2% over a region that 
experiences an atom-density variation of 35%. 

III. VORTEX CORES 
A. Core size 

The other defining length scale of the vortex lattice is 
the core radius. Here we study the core radius in the 
Thomas-Fermi regime (as opposed to the lowest Landau 
level regime, described later) where it should scale with 
the healing length. A theoretical value for the vortex core 
radius was generated by performing a numerical simu- 
lation for a 3D BEC containing an isolated vortex and 
comparing the fitted radius of this vortex to the corre- 
sponding healing length. Fitting the simulation in the 
same manner that we later treat the experimental data 
(described in Sec. I B) we obtain an expression for the 
core radius of 

r v = 1.94 x £ , (4) 

with healing length £=(8-7ma sc ) -1 / 2 , where a sc is the 
scattering length and n is the density- weighted atom den- 
sity. For the data presented, n is determined from the 
in-trap image before expansion. 

Core size measurements and fractional core area (dis- 
cussed in the next subsection) measurements require con- 
siderable attention to detail. In pursuing these measure- 
ments, we find that nearly everything — from focusing 
issues, to lensing due to off resonant imaging light, to 
even imperfect atom transfer into the anti-trapped state 
before expansion — can lead to an overestimation of the 
vortex core size. By far the biggest potential systematic 
error in our system is axial expansion, which, as noted in 
Sec. I B, requires careful attention. 

A range of core sizes is achieved by varying the initial 
number of atoms loaded into the magnetic trap prior to 
evaporation. To avoid the core size saturation effect, due 
to high condensate rotation , we consider only clouds 
with Tlll > 10, where Tlll = ^/{2Ml) is the LLL 
parameter and fj, is the chemical potential. This ratio 
of chemical potential to rotation approaches unity as we 
enter the LLL regime, while at values of 10 or greater 
we should be firmly in the Thomas-Fermi regime. In 
practice this requires only that we keep the condensate 
rotation rate low. Core size is measured by fitting the 
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FIG. 3: Comparison of measured core radii with the Thomas- 
Fermi prediction (Eq. @) represented by the solid line. Black 
squares are the core size in expansion scaled by the radial 
expansion of the condensate so that they correspond to the 
in-trap values. Data shows reasonable agreement with theory. 
The fact that the measured core size is consistently larger is 
likely due to the fact that nearly all our imaging systematics 
lead to an overestimation of the vortex core size. 



expanded image with a Thomas-Fermi profile and each 
core with a 2D Gaussian. For Fig. [3| the measured core 
radius in expansion is scaled back to the corresponding 
in-trap value using the radial expansion factor discussed 
in Sec. I B. To reduce scatter we consider only vortices 
located less than half a condensate radius out from the 
center. Additionally we find that some vortices appear 
to have some excitation or bending which leads to a poor 
fit. To filter these out we consider only vortices that have 
a contrast greater than 0.6. Here contrast is defined, 
with respect to the integrated (along the line of sight) 
condensate profile, as the peak of the "missing" column 
density at the vortex position divided by the smoothed 
Thomas-Fermi profile at the same position. On average, 
about 30% of visible vortices fulfill all the criteria for 
being used in the core size measurement, although this 
number can vary wildly (10%-100%) depending on vortex 
number and position. 

From Fig. |3| we can see that the data and the Thomas- 
Fermi theory agree reasonably well. The data do seem 
to be slightly above the theory value, but we are hesitant 
to make too much of this because, as noted before, there 
are many systematic errors that tend to bias the data 
toward larger core size. Measurement is easier and the 
agreement better, on the low density, large-core side of 
the graph. 

At an early stage in this work, we speculated that 
the mean-field Gross-Pitaevskii equation might not give 
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a good quantitative description of vortex core size be- 
cause the core size is particularly sensitive to the healing 
length £. At our highest densities, while the gas is nom- 
inally dilute (na? sc < 10 -5 , where a sc is the interatomic 
scattering length), the mean interatomic distance rt -1 / 3 
is only a factor of 1.5 less than £. Our data, however, 
are ambiguous with respect to the accuracy of the Gross- 
Pitaevskii equation for predicting vortex core sizes at our 
highest densities. The roughly 25% discrepancy between 
our measurements and the mean-field prediction shown 
at the smaller radii in Figgis comparable to possible sys- 
tematic errors in our measurements of the smaller cores 
that exist at higher densities. In retrospect, our exper- 
imental design is such that we may be unable to see a 
mean- field failure even if one were to exist. During the 
radial expansion, the density drops. Thus the accuracy 
of the mean-field approximation is likely to improve sig- 
nificantly during the expansion. Our anti-trapped expan- 
sion, while more rapid than a conventional ballistic ex- 
pansion, is still slow compared to the rate at which a vor- 
tex can adiabatically relax its radius |2fij (approximately 
fJ-/h). Any non-mean-field corrections to the vortex core 
size will likely relax away before the cores have ex pan ded 
to be large enough for us to reliably image them 36] . 



densate enters the LLL regime. An alternate treatment 
from Baym and Pethick |25j predicts that A saturates at 
0.225 as the vortices go from a Thomas-Fermi profile to 
the profile of a LLL wave function. Our numerical data 
for A, together with experimental points, are plotted in 
Fig. E|(a,b). For T^ LL < 0.1, the data agree reasonably 
well with the Thomas-Fermi result. For larger V^ LL , the 
data clearly show a smooth transition to the LLL regime, 
and a saturation of A at the LLL limit. 

The experimental data in Fig. Efb) tend to lie above 
the numerical data. This is likely due to our tendency 
to overestimate the core size, as discussed previously in 
Sec. Ill A. Figure 0Jc) demonstrates the dangers of ax- 
ial expansion in this measurement. For the data pre- 
sented the condensate undergoes a factor of 2-3 in axial 
expansion, and we see a corresponding increase in A. 
This clearly illustrates the importance in suppressing ax- 
ial expansion for these measurements. It is interesting to 
note that with our rapid axial expansion the fractional 
core area can overshoot the LLL limit, which in principle 
should still be valid in the limit of adiabatic expansion. 

C. Core density profile 



B. Fractional core area 

As noted in the introduction, the fractional conden- 
sate area occupied by the vortex cores is also a quan tity 
that has been of much theoretical interest [lfi. 1 2.4 l27j 
and has been experimentally studied previously [29j. It 
is argued by Baym and Pethick [2^ that the fractional 
core area reaches a limiting value as one enters the LLL 
regime. The corollary of this argument is that fractional 
core area is a reasonable way to monitor the transition to 
the lowest Landau level regime. We examine this transi- 
tion with numerical work, which we can to push further 
into this regime than we can achieve experimentally. Ad- 
ditionally we examine some of the systematic errors that 
can affect the experimental data. To this end numeri- 
cal calculations were performed as previously described, 
for 3 x 10 6 ,5 x 10 5 ,and 1 x 10 5 atoms, and for rotations 
ranging from il/uj p =0.15 to 0.998. For the experimental 
data, actual condensates were generated over a similar 
range with Q/uj p =0.15 to 0.98 and N=4- 50 x 10 5 . The 
numerical data as well as the experimental data are fit in 
the same manner as described in Sec. I B. 

We define the fractional area A occupied by the vor- 
tices to be A — n v irr 2 , where n v is the areal density of 
vortices and r 2 is defined in Eq. 0] Ignoring density in- 
homogeneity effects, in the limit of many vortices, the 
expected vortex density n v is mfl / (ttK) . The resulting 
prediction for A can be expressed as A = 1.34 x T^ LL . 
This value exceeds unity for Tlll < 1-34, which has led 
to the prediction that vortices should merge as the con- 



We can also observe the transition to the LLL regime 
in the numerical data by examining the shape of the 
condensate vortex cores. In the Thomas-Fermi regime, 
the vortex-core density profile is well described by the 
form n(r) = (rl \/l( 2 + r 2 ) 2 |2oj where r is measured 
from the vortex center. Alternatively in the LLL regime, 
the core is no longer dictated by the interactions but 
rather by kinetic energy considerations. In this case, 
within the Wigner-Seitz unit cell, the vortex is thought 
[25| to have a simple oscillator p-state structure n(r) = 
((Cr/6) • exp[-r 2 /2l 2 ]) 2 for < r < I, where C is a nor- 
malization constant, and I is the radius of the Wigner- 
Seitz unit cell and is related to the nearest-neighbor lat- 
tice spacing b by I = (v / 3/27r) 1 / 2 fe. Figure|S]is a compar- 
ison of the central vortex, in three numerically generated 
condensates, to both Thomas-Fermi and LLL predicted 
core shapes. The simulation for Fig. [5Ja) was performed 
for 3 x 10 6 atoms and Q,/lo p = 0.15 and is well inside 
the Thomas-Fermi regime (Tlll = 117). Here the den- 
sity profile of the numerical data (solid line) seems to 
fit quite well to the Thomas-Fermi vortex form (dotted 
line), but the LLL form is a poor description of the vortex 
core (dashed line). The simulation for Fig.|SJb) was per- 
formed for 5 x 10 5 atoms, ^l/u p — 0.95 and Tlll = 3.6. 
One can see from Fig. 0] that this is in the transition 
region. Not surprisingly both vortex forms fit about 
equally well. In Figs. Efb) and Fig. OJc), the vertical 
line represents the edge of the Wigner-Seitz unit cell at 
r = I. The simulation for Fig. E{c) was performed for 
1 x 10 5 atoms, fl/uj p = 0.998 and T L ll = 0.72. One can 
see that LLL is a much better description of the vortex. 
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FIG. 4: Fractional condensate area occupied by vortex cores 
(A) as a function of T~[ L l = 2Hfl/fX, the inverse lowest Lan- 
dau level parameter. Plot (a) shows a smooth transition in 
the numerical data from the Thomas-Fermi limit where A is 
linear in T~[ LL to the LLL limit where A saturates. Here the 
Thomas-Fermi theory is represented by the dashed line and 
the LLL limit by the dotted line. Plot (b) is a comparison 
of the numerical data to experimental data. Plot (c) demon- 
strates the effect on the experimental measurement of allow- 
ing the condensate to expand axially during the expansion 
process. 
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FIG. 5: Numerically generated vortex core density profiles 
approaching the lowest Landau level regime. Solid lines rep- 
resent the numerical result for (a) Tlll = 117, {V)Tlll = 
3.6, (c)Tlll ~ 0.72. The dashed line is the expected profile 
for a LLL wave function 25] given the condensate rotation. 
The dotted line is the expected vortex form in the Thomas- 
Fermi limit 23 given the condensate density. The vertical 
lines in figure (b) and (c) designate the edge of the Wigner- 
Seitz unit cell. As Tlll decreases, one can see a clear transi- 
tion from the interaction-dominated Thomas-Fermi regime to 
a LLL function where kinetic energy concerns and the vortex 
core spacing dictate the shape and size of the vortex. 
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On a separate but interesting note, even this far into 
the LLL regime, our numerical solution of the GP equa- 
tion shows that the radial profile of the overall smoothed 
condensate fits much better to a parabola than to the 
Gaussian that was originally predicted 0. The reason 
the Gaussian-density-profile prediction fails to pan out 
can be extrapolated from data presented earlier in this 
paper. The density-profile prediction for the radial pro- 
file in the LLL arose from an elegant argument that was 
based on an assumption that the vortex nodes were on 
a perfect triangular lattice. As was originally pointed 
out to us by A. H. MacDonald j3?| and has been the 
subject of two recent theoretical works [Isl flf|. a slight 
radially dependant perturbation in the areal density of 
vortices is enough to convert a Gaussian density distribu- 
tion into an inverted parabola. The analytic description 
of this perturbation in (calculated in the LLL) bears 
a striking resemblance to the one measured in Sec. II in 
the Thomas-Fermi regime and also to the analytic form 
^3 calculated in the Thomas-Fermi regime. The sur- 
prising result of this perturbation in the areal density 
of vortices is that one of the most striking features of 
the Thomas-Fermi regime, the parabolic Thomas-Fermi 
density profile, still exists in the LLL regime where the 
condensate kinetic energy is clearly non-negligible com- 
pared to interaction energy. 



IV. ROTATIONAL SUPPRESSION OF THE 
QUANTUM DEGENERACY TEMPERATURE 

So far in this article, condensate rotation has been con- 
sidered in terms of its discrete vortex elements and their 
dynamics. A more global effect of condensate rotation 
is to suppress the quantum degeneracy temperature, T c , 
for a fixed number of trapped particles, N |2^|. While 
superficially similar to magnetic- field suppression of T c in 
superconductors, this effect is distinct in that it is not a 
many-body effect but rather a suppression of atom den- 
sity due to a rotationally weakened trapping potential. 

One can define a static critical temperature, that 
applies in the case of no rotation. For a non-interacting 
gas in a harmonic potential, = 0.94?iw^ o A 1 ' /3 , 

where N is the number of atoms in the sample and 
ujho = (w^Wz) 1 / 3 is the average trap frequency. Rota- 
tional suppression of T c may be accounted for by the cen- 
trifugal weakening of the radial harmonic confinement, 
1 — fi 2 /wp, which leads to a reduction in the 

thermal gas density n t h relative to the non-rotating case. 
As a result, for fixed N, a lower temperature must be 
reached before the phase space density (~ n t h/T 3 ^ 2 ) is 
sufficiently high to bring about the BEC transition. The 
expression for T c as a function of the rigid-body rotation 
rate f2 of the gas sample is 
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FIG. 6: (Color online) Condensate fraction versus tempera- 
ture for various sample rotation rates. Condensate number, 
Nbec, and total number, N, are obtained from fits to cloud 
images. The temperature T is extracted from the vertical 
width of the thermal cloud while the centrifugal distortion of 
the thermal cloud's aspect ratio yields its rotation Q n- The 
temperature is scaled by the static critical temperature Ti°^ 
for an ideal gas. The data has been grouped according to 
three different ranges of Qn- The solid line is the theoretical 
dependance expected for a static ideal gas. 



In Fig. El experimental data have been used to plot con- 
densate fraction versus temperature for three different 
ranges of sample rotation rate. All quantities are ob- 
tained from fits to the non-destructive, in-trap images of 
the trapped gas clouds. As usual, the rotation rate has 
been assessed from the changing aspect ratio of the ther- 
mal cloud according to Eq. ^ It is qualitatively clear 
from Fig. [S] that a lower temperature is required at 
higher rotation rates to reach a given condensate fraction. 
For each data point, the temperature, T, and conden- 
sate fraction, Nbec/N (where Nbec is the number in 
the condensate) can be used together to infer a value for 
the critical temperature of the sample using the formula 
T c = [1 - Nbec/N]~ 1/3 T. To remove shot-to-shot vari- 
ation in T c due to atom-number fluctuations, the value 
of the inferred transition temperature can be scaled by 
the static value calculated from the measured atom 
number. The scaled value of T c is plotted against rota- 
tion rate in Fig.[7|for three ranges of temperature: "hot," 
"medium" and "cold," corresponding to three different 
ranges of condensate fraction as given by the legend of 
Fig.0 

The centrifugal suppression of T c is perhaps less inter- 
esting than suppression due to many-body interactions. 
Nevertheless, from a technical stand point, centrifugal ef- 
fects are an important consideration in the evaporative 
cooling and spin-up process. 
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FIG. 7: (Color online) Inferred critical temperature T c , scaled 
by the non-rotating expectation Tc f° r an ideal gas, as a 
function of thermal gas rotation For each data point, T c 
is inferred from the measured condensate fraction and tem- 
perature of a sample. The total number of atoms is used to 
obtain Tj°\ The data has been grouped according to three 
different condensate fractions as shown. A set of static points, 
where the thermal cloud is not stirred before evaporation, are 
deliberately plotted at zero rotation. Otherwise, the rotation 
rate fijv is obtained from the thermal clouds aspect ratio. 
The rotating data has been cropped at a minimum threshold 
0.2 < Sliv/cjp to avoid imaginary rotation values arising from 
noise in near-static aspect ratios. The dotted line is the the- 
oretical expectation according to Eq. with Tj ' as for an 
ideal gas. The solid line is a fit of the data to Eq. with 
an arbitrary overall scaling of the vertical axis. The fit result 
is equivalent to assuming an effective that is 87 ± 1% of 
the ideal gas value. This 13% shift is consistent with what 
one would expect from atom-atom interactions |3lt . 



V. CORE CONTRAST AND CONDENSATE 
TEMPERATURE 

Since the very first observations of dilute-gas BEC, the 
temperature of the sample has been determined by imag- 
ing the "skirt" of thermal atoms that extends beyond the 
radius of the condensate. In practice, it is difficult to ex- 
tend this measurement below about T/T c — 0.4, except 
in very special cases (for instance when a Feshbach res- 
onance is used to set the scattering length to zero). Fol- 
low temperatures, the density of thermal atoms becomes 
so low that they are difficult to image. Moreover, when 
the temperature becomes lower than or comparable to 
the chemical potential of the self-interacting condensate, 
the spatial extent of the thermal cloud is no longer ap- 
preciably larger than the condensate itself. 

It was suggested that vortex cores might serve as 
"thermal- atom concentration pits" , in order to enhance 



thermometry at low temperatures. In a simple Hartree- 
Fock (HF) picture of the interaction between thermal 
atoms and the condensate, the condensate density rep- 
resents a repulsive interaction potential to the thermal 
atoms. Along the nodal line of a vortex core, the con- 
densate density and presumably its repulsive interaction 
potential vanish. Thus, the thermal atoms would ex- 
perience the lowest combined interaction and magnetic 
potential within the cores of vortices. As a result their 
density would be highest there. Additionally, images of 
thermal atoms in the vortex core could be taken against 
a vanishing background condensate density. Moving be- 
yond the HF approximation, one finds a more compli- 
cated picture. The Bogoliubov spectrum of very long 
wave-length thermal phonons extends all the way down 
to the chemical potential. One should contrast this en- 
ergy with the energy of a thermal atom confined to a 
vortex core. Perhaps the atom experiences no interaction 
energy. However, the kinetic energy cost of bending its 
wave function to fit inside a core with a radius of the or- 
der of healing length must, by definition, be comparable 
to the chemical potential. In the limit of very elongated 
vortex cores, there can be very low-energy, core-bending 
modes |5j, |39( . Thermal excitations of these modes would 
manifest as a temperature-dependent contrast ratio. We 
expect this effect is unlikely to be important in the rela- 
tively flattened geometry of our highly rotating conden- 
sates. In any case, without more rigorous analysis, it is 
not easy to predict how the contrast ratio of our vortices 
should vary with temperature, but we nonetheless set out 
to do a preliminary study of the effect. 

We vary the final condensate temperature by chang- 
ing our rf-evaporation end point. This produces a cloud 
with temperatures between 5 — 50 nK or T/T c between 
1 and less than 0.4. Here T c is calculated from the 
trap frequencies and a measurement of total atom num- 
ber using the formula T c — 0.94/i(I)^ o iV 1 / 3 , where u>ho 
has been adjusted for rotation according to the equation 
Cuho = ^/io(l — ^ 2 / w p) 1 ^ 3 - When possible, T is extracted 
from a two-component fit to the in-trap image. Because 
our rotation rate and temperature are linked through 
the ID evaporative process, it is unavoidable that also 
varies during the data set. 

To measure core contrast, we expand the cloud using 
the usual expansion procedure. The atom cloud is ex- 
panded radially by a factor of 13 to ensure that the cores 
are large compared to our imaging resolution. However, 
because we no longer care about the precise core size 
we do not suppress the axial expansion. Additionally, 
the axial expansion actually reduces background fluctua- 
tions in the measured core contrast. With a factor of two 
axial expansion, cores become much rounder and clearer 
as shown in Fig.^c). These changes allow us to achieve 
a higher core contrast and quieter signal than we can 
without expansion. 

The term core brightness (1-contrast ratio) will be our 
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metric for this experiment. We define core brightness 
(B) as ri2D (core) /ri2D (cloud) , where riiD(core) is the ob- 
served atom density, integrated along the line of sight, 
at the core center, and n 2 D(cloud) is the projected in- 
tegrated atom density at the same point, based on a 
smoothed fit to the overall atom cloud. To determine 
n2D(cloud), we fit the condensate image to a Thomas- 
Fermi profile and the surrounding thermal atoms to a 
Bose distribution. We find n2u(core) by fitting each vor- 
tex with a Gaussian to determine its center and then 
averaging five pixels around the center point to deter- 
mine the integrated density. Brightness is calculated for 
each vortex and then averaged with other vortices in the 
cloud. To suppress noise from low signal, vortices further 
than 0.4 R p from the condensate center are disregarded 
for this measurement. The ri2D(core) term necessarily 
contains signal from the surrounding thermal atoms be- 
cause the vortices do not penetrate the thermal compo- 
nent. Thus, one expects to see a steady decrease in B 
with decreasing temperature, as atoms not necessarily in 
the vortex core, but still in the integrated line of sight, 
disappear. One would hope that B continues to decrease 
even for T/T c below 0.4 for this analysis to be a viable 
means of extending condensate thermometry. 

We are in the awkward position of comparing 
our core contrast measurement to a temperature 
measurement that, as previously described, is ex- 
pected to fail at low temperatures. To monitor 
this failure, we calculate a simplistic core brightness 
(B simple) found by comparing the fitted in-trap con- 
densate and thermal cloud profiles. Here B S i mp i e = 
fi2D (thermal) / (fi2D (condensate) +fi2D (thermal)) where 
fi2D(condensate) and n2D (thermal) are the smoothed 
condensate and thermal cloud profiles integrated along 
the z-axis and averaged over a region of radius less than 
0.4 R p from the condensate center. The term B S i mp i e 
can be thought of as the core brightness one would ex- 
pect based on the undoubtedly false assumption that the 
condensate and thermal atoms do not interact. It is in- 
teresting to compare B to B S i mp i e since this same dubi- 
ous assumption is implicit in the standard thermometry 
technique of fitting the thermal "skirt" . 

In Fig. HO B and B simp [ e are plotted versus the final 
evaporative cut. For our experiment, the thermal cloud 
can be reliably fit for T/T c > 0.6 and less reliably fit 
for T/T c > 0.4. In both these regions T/T c decreases 
continuously with lower final evaporative cut. It is as- 
sumed that for T/T c just below 0.4, this trend contin- 
ues. For reference, three values of T/T c (measured from 
the thermal "skirt") are included in the plot. One can 
see that B does steadily decrease with lower temperature 
for T/T c > 0.4. It is interesting to note that B S i mp i e 
closely tracks B at the higher temperatures and then 
diverges from B as the cloud gets colder. Presumably, 
this divergence occurs because thermal atoms are pushed 
away from the condensate center as interactions between 
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FIG. 8: Measured core brightness as a function of final rf evap- 
orative cut. Within our ability to measure, T/T c decreases 
continuously with the rf frequency. For the black squares 
brightness (B) is defined as the 2D atom density at the vortex 
core divided by the 2D atom density of the overall smoothed 
condensate plus thermal cloud profile at the same point. For 
the open triangles a simplistic brightness (B s i mp i e ) is calcu- 
lated from the ratio of the 2D atom density of the thermal 
cloud to the 2D atom density of the overall smoothed con- 
densate and thermal cloud profile. At high temperatures B 
and B S i mp i e exhibit a clear dependance on the final rf cut. At 
lower temperatures it is encouraging that as £> s ; mp ; e , begins 
to fail B is still continuing a smooth trend downward. Dis- 
appointingly at very low temperatures, B plateaus at about 
0.14. 



the condensate and the thermal cloud become important. 
The fact that B S i mp \ e diverges upwards is likely due to 
the tendency of our fitting technique to overestimate the 
thermal cloud density at high condensate fractions. The 
failure of B S i mp i e at low temperatures also throws into 
suspicion the quoted T/T c since they are determined from 
the same two-component fit. 

In contrast, as B S i mp i e begins to fail, B continues its 
previous smooth downward trend. It is also interesting 
to note that at an rf of 2.35 MHz, we see a B of 0.13-0.15, 
which is not that far off from the work of Virtanen et al. 
[22]| who predict that atoms trapped in the core would 
lead to a B of 0.1 at a T/T c of 0.39. Unfortunately, our 
efforts to observe a B of less than 0.125 have failed so far, 
as can be seen from the data points at 2.3 MHz in Fig. El 
This limit impedes our ability to measure temperatures 
colder than 0.4 T/T c . Currently, it is unclear what the 
source of this limit is. Perhaps the same imaging system- 
atics that make our vortex radius unreliable at the 10% 
level are also preventing us from seeing a core brightness 
level less than 0.13 or a very slight tilt of the vortices 
may occur during expansion. 



12 



As a caveat to the previous discussion, the same limita- 
tions that inhibit condensate thermometry below T/T c of 
0.4 will also reduce the efficacy of evaporative cooling in 
the same regime. Additionally, the already inefficient ID 
nature of our evaporation would exacerbate such a cool- 
ing problem. Perhaps the simplest explanation for the 
failure of B to decrease with very deep rf cuts is that the 
condensate fraction is no longer increasing. One could 
imagine the our measured B is faithfully following the 
temperature we achieve. 

In summary the conclusions of our preliminary attempt 
to extend thermometry with core brightness are encour- 
aging but ambiguous. New ideas are needed before we 
can make further progress. 

The work presented in this paper was funded by NSF 
and NIST. We acknowledge fruitful discussions with Dan 
Sheehy and Leo Radzihovsky. Numerical simulations 
were performed with the JILA Keck cluster. 
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